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I.  INTRODUCTION 


The  purpose  of  this  research  is  to  make  a  detailed  investigation 
of  finite-amplitude  standing  waves  in  air  at  ambient  temperatures  in 
a  rigid  walled  cavity,  for  comparison  with  the  theory  oi  Coppens  and 
Sanders  [l]  modified  to  include  observed  energy  losses.  This  in¬ 
vestigation  is  a  continuation  of  the  experimental  works  of  Beech  [2] 
and  Winn  C 33 ,  with  an  extension  to  a  rectangular  cavity  vice  a  long 
round  tube. 

A  rectangular  cavity  was  chosen  as  it  was  a  step  toward  the 
study  of  two  and  three  dimensional  waves,  although  this  investi¬ 
gation  concerns  only  one  dimensional  standing  waves.  In  addition, 
it  was  thought  that  a  rectangular  geometry  would  produce  over¬ 
tones  which  would  be  less  harmonic  than  those  for  the  more  sym¬ 
metric  geometry  of  the  round  tube.  The  dimensions  of  the  cavity 
were  chosen  in  an  effort  to  space  out  the  system  resonances  to 
avoid  degeneracy  or  near -degeneracy .  This  was  accomplished  with 
some  success,  and  thus  allowed  for  the  accurate  location  and  anal¬ 
ysis  of  the  first  few  overtones  of  each  of  the  two  modes  studied. 


II.  BACKGROUND  AND  THEORY 


A  finite  amplitude  wave  is  a  solution  of  the  non-linear  wave 
equation  that  describes  the  real-life  environment  that  confronts 
large  amplitude  acoustic  waves.  This  finite-amplitude  wave  re¬ 
quires  second  order  terms,  or  higher,  to  describe  the  waveform. 

It  has  been  known  for  many  years  that  hydrodynamic  equations 
describing  traveling  waves  in  a  non-dissipative  medium  predict  a 
change  in  the  waveform.  Kirchoff  [4]  and  Lamb  [5]  were  among  the 
early  investigators  in  this  area.  In  more  recent  years  there  have 
been  many  theories  and  methods  of  solution  proposed  to  describe 
this  distorted  traveling  waveform  in  a  non-dissipative  medium. 
Amoung  these  are  the  works  of  Fay  [6],  and  Blackstock  [7]  who  pro¬ 
duced  an  extensive  theoretical  treatment  of  low  amplitude  non¬ 
linear  waves,  and  made  use  of  Fourier  and  power  series 
approximations . 

When  the  problem  was  extended  to  dissipative  mediums,  it  was 
shown  very  early  [8]  that  waves  in  absorptive  media  exhibit  finite- 
amplitude  characteristics  only  at  high  strengths,  different  from 
waves  in  a  non-dissipative  media.  More  recently,  an  important  con¬ 
tribution  to  the  understanding  of  finite  amplitude  waves  in  a  dis¬ 
sipative  medium  was  presented  by  Fay  [9].  He,  similar  to  a  later 
work  by  Blackstock  [lO],  formulated  a  Fourier  series  solution  to 
describe  the  finite-amplitude  waveform  for  a  plane  wave  of  infinite 
extent.  More  recent  theoretical  appioaches  to  finite-amplitude 
effects  have  been  of  two  basic  types:  an  approximation  solution, 
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such  as  that  proposed  by  Fox  and  Wallace  Cll],  and  a  perturbation 
solution  like  that  of  Fay  C9 ] .  The  approximation  solution  is  based 
on  assumed  attenuation  mechanisms  and  manner  of  overtone  population. 
The  perturbation  solution  makes  use  of  perturbation  methods  and 
Fourier  series  representation  of  the  waveform  C 10, 12,13].  A  per¬ 
turbation  analysis  which  was  extended  to  the  sixth  order  by  Keck 
and  Beyer  C 133  is  of  particular  importance  to  this  investigation, 
as  will  be  evident  later .  The  solution  obtained  by  Keck  and  Beyer 
is  valid  near  the  source,  unlike  previous  perturbation  solutions 
of  finite-amplitude  traveling  waves  in  dissipative  media. 

Although  most  investigations  have  been  confined  to  traveling 
waves,  there  has  been  some  investigations  of  standing  waves.  For 
example,  Keller  [l4]  extracted  solutions  for  standing  waves  in  a 
closed  tube.  In  that  investigation,  Keller  assumed  a  non-dissi- 
pative  medium  which  resulted  in  infinite  amplitude  waves  at  re¬ 
sonance,  yet  still  obtained  useful  results  for  frequencies  close 
to  those  of  resonance.  Coppens  and  Sanders  Cl]  developed  an  ex¬ 
tension  of  the  Keck-Beyer  perturbation  approach,  that  include  wail 
losses  as  predicted  by  the  Rayleigh-Kir choff  theory  [15].  This 
approach  was  carried  to  the  tenth  order  by  Ruff  [l6],  at  which 
point  it  was  found  that  the  mathematics  become  divergent. 

Coppens  [l7]  later  expanded  the  theory  and  used  a  Fourier  decom¬ 
position  technique,  which  both  Beech  C 3*3  and  Winn  12]  experimentally 
investigated.  It  was  found  that  the  tubes  that  they  investigated 
did  not  behave  in  the  manner  predicted  by  Rayleigh-Kir choff  loss 
mechanisms.  This  investigation  follows  up  the  observation  that 
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the  tubes  exhibited  non-Ray leigh-Kir choff  behavior,  by  testing  the 
theory  using  empirically  determined  losses. 

In  the  theory  of  Coppens  and  Sanders,  a  one  dimensional,  non¬ 
linear  wave  equation  with  dissipative  terms  which  correspond  to 
those  encountered  in  the  discription  of  plane  standing  waves  in  a 
rigid-walled  cavity,  is  developed.  The  wave  equation  in  this 
theory  can  be  Fourier  decomposed  to  take  the  form 


k2- 


n=l\bx2  C2  bt2 
o 


+  D 


=  b  lV 


n  I  bx 


bx' 


2  l  bx 


(1) 


where 


D  =  6 
n  1 


_1_ 

r. 


3/2  v  2V  %  .  2 
CO.^nJ  bx  bt  n  bx 


and  contains  a  dispersive  term  and  a  dissipative  terra,  both  based 
on  the  Rayleigh-Kir choff  theory  of  wall  losses  [15],  and 

h  ■ 

§n  ~  particle  displacement  due  to  harmonic 

5  =  s?n 

n  n 

b  =  (1  +  Cp/Cv)/2 

Cp  =  specific  hear  of  medium  at  constant  pressure 
Cy  =  specific  heat  of  medium  at  constant  volume 


CO. 


=  frequency  (angulc.r)  of  fundamental  resonance 


CO. 


Q1  ~  CO. 


lr 


1+ 


CO. 


1- 


and  C0^__  are  the  half  power  frequencies  above  and  below  the 
fundamental . 

In  the  description  of  waves  in  cavities,  Eq.  1  can  be  refor¬ 
mulated  [l8j  to  introduce  the  actual  cavity  resonances  and 
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-  peak  pressure  of  fundamental  component  of  wave 
Rn  =  Four  iei'  coefficient  of  n  harmonic  component,  normalized 
such  thst  R^  =  1 

and 

0n  =  phase  angle  of  n^*1  harmonic  component,  where  0^=0. 
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This  equation  can  be  manipulated  into  a  general,  time-independent 
form  which,  when  evaluated  at  the  rigid  boundaries  of  the  cavity, 
can  be  furthe*  simplified  to  produce  the  following  form  amenable  to 
computer  calculation: 


H  R  COS  (0  -0  ) 
n  n  '  n  n' 


MbQx 


R  .R  COS  (0  .  + 
3  n-J  3 


0  •) 

n-J 


R  ,  .R  COS  (0  . 

n+j  j  '  n+j 


and  a  similar  form  with  "sin"  in  place  of 


where 


"aw", 


and 


2ACU  Q,/C0  . 

tan  9  =  ■ 

n  0,/Q,, 


As  can  be  seen  from  the  above  equations,  it  is  important  to  know 
the  Q(N) ’s  and  W  's  of  the  system  for  the  input  of  the  theoretical 
predictions  calculated  by  the  computer  .  It  is  convenient  to  de¬ 
fine  a  quantity,  E(n),  that  indicates  the  position  of  CO  relative 
to  the  classical  harmonic  frequencies,  n(0^r .  E(n)  is  defined  by 

E(n)  =  (<0nr/(nt0lr ) )  -  1. 

The  values  of  Q(n)  and  E{n)  are  directly  determined  from  the 
infinitesimal  amplitude  analysis  of  the  cavity  and  are  the  input 
quantities  for  the  computer  calculations. 
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III.  EXPERIMENTAL  CONSIDERATIONS 


A.  APPARATUS 

A  block  diagram  of  the  experimental  apparatus  is  shown  in 
Figure  1.  Finite-amplitude  waves  were  generated  in  air  at  ambient 
temperatures  in  a  rectangular  cavity  with  interior  dimen', '.ons  of 
6.97  cm.  by  20.96  cm.  by  30*28  cm.  The  walls  of  the  ca.ity  were 
constructed  of  0.985  in.  thick  aluminum  plates,  milled  and  screwed 
together  to  provide  t-ght,  right  angle  joints.  These  joints  sealed 
by  applying  thin  layers  of  silicon  grease  to  adjoining  surfaces 
before  screwing  them  together.  Acoustic  waves  were  generated  in 
the  cavity  by  exciting  the  air  with  a  plane-faced  piston.  The 
piston  face  gained  access  to  the  interior  of  the  cavity  by  way  of 
an  orifice  bored  in  one  of  the  cavity  walls.  The  piston  was  fitted 
with  two  O-rings  and  the  piston  orifice  was  bored  such  that  the 
piston  ivith  0-rings  would  make  an  airtight  fit.  It  was  found  that 
the  piston  motion  was  restricted  by  the  second  0-ring;  therefore 
the  piston  was  inserted  into  the  orifice  only  as  far  as  the  first 
0-riug.  Lubrication  and  further  sealing  of  the  piston-orifice 
interface  was  provided  by  the  same  high  vacuum  silicon  grease  used 
to  seal  the  wall  joints.  The  piston  was  attached  to  an  M-B  Elec¬ 
tronics  Model  EA  1500  Exciter.  This  exciter  was  powered  by  two 
M-B  Electronics  Model  2120MB  Power  Amplifiers  operating  in  parallel. 
The  maximum  attainable  acceleration  of  the  piston  with  this  equip¬ 
ment  was  about  65  g*s  in  free  air  at  100  Hz.  The  frequency  and 
power  level  of  the  piston  were  controlled  by  a  General  Radio  Il6l-A 


11 


tre  1.  Apparatus 


Coherent  Decade  Frequency  Synthesizer.  Frequencies  can  be  selected 
to  within  less  than  0,01  Hz  with  this  synthesizer.  The  piston  mo¬ 
tion  was  monitored  by  measuring  the  output  of  an  Endeveco  Model  2215 
Accelerometer  mounted  within  the  piston.  Distortion  in  the  piston 
motion,  as  determined  by  the  accelerometer,  set  the  lower  limit  to 
which  reliable  measurements  of  the  non-linear  distortion  in  the 
pressure  waveform  could  be  made.  Measurement  of  the  pressure  wave¬ 
form  in  the  cavity  was  made  in  the  far  corner  from  the  piston 
source.  Refer  to  Figure  2.  The  sound  pressure  levels  and  harmonic 
components  of  the  pressure  wave  in  the  cavity  were  monitored  by  a 
1/4  inch  diameter  Bruel  and  Kjaer  type  4136  Condenser  Microphone. 
This  microphone  was  mounted  such  that  Ihe  active  monitoring  element 
of  the  microphone  was  flush  with  the  interior  of  the  cavity.  The 
microphone  was  fitted  with  an  0-ring  and  seated  in  its  orifice  with 
the  same  high-vacuum  silicon  grease.  The  output  signal  of  the 
microphone  was  amplified  by  a  Bruel  and  Kjaer  type  2801  Power 
Supply.  This  power  supply  with  the  microphone  comprized  the  mi¬ 
crophone  system.  The  microphone  system  was  monitored  by  a  Hewlett- 
Packard  Model  400D  Vacuum  Tube  Voltmeter.  Another  vacuum  tube 
voltmeter  of  the  same  type  was  used  to  monitor  the  output  of  the 
accelerometer.  A  Hewlett-Packard  Model  302A  Wave  Analyzer  was 
used  to  measure  harmonic  distortion  of  both  the  accelerometer  and 
microphone  output  signals.  A  Hewlett-Packard  Model  521A  Electronic 
Counter  was  used  to  monitor  the  filtered  output  of  the  wave  ana¬ 
lyzer  to  assure  that  the  proper  frequency  was  being  analyzed. 
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Microphone  Orifice 


B.  EQUIPMENT  CALIBRATION 

Absolute  calibration  of  the  microphone  system  response  was 
required  to  determine  strength  parameters.  The  microphone  system 
included  the  microphone,  its  preamplifier,  and  the  vacuum  tube 
voltmeter  connected  to  the  output  of  the  preamplifier.  The  micro¬ 
phone  system  response  is  defined  by 

S  =  V  /P 
m  m 

where  Vm  is  the  RMS  voltage  from  the  preamplifier  as  indicated  on 
the  VTVM  and  P  is  the  RMS  pressure  amplitude  at  the  microphone 
face.  Calibration  of  the  microphone  system  was  accomplished  with 
the  use  of  a  Breul  and  Kjaer  Model  4220  Pistonphone  which  pioduces 
a  124  i  0.2  re  volt/p-bar  acoustic  signal  at  250  Hz.  The  micro¬ 
phone  system  sensitivity  was  found  to  be 

sm  =  U*46  +0,03)  10"3  volt/(Nt/m2) , 

which  is  in  close  agreement  with 

sra  =  1*50  +  0.03  x  10~3  Volt/ ( Nt/m2 ) 

and 

Sm  =  1.52  +0.03  x  10~3  Volt/(Nt/m2) 

as  determined  by  Beech  [2]  and  Winn  [3]  respectively.  The  micro¬ 
phone  system  had  a  flat  frequency  response  to  within  +  0.5  db 
from  80  Hz.  to  20  kHz.,  as  determined  from  the  microphone  response 
curve  supplied  by  the  manufacturer.  With  the  microphone  system 
calibrated,  the  wave  analyzers  were  calibrated  with  respect  to  the 
microphone  system  VTVM.  This  was  accomplished  by  removing  the 
microphone  and  its  preamplifier  from  the  VTVM,  and  substituting 
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a  voltage  signal  from  the  frequency  synthesizer  to  both  the  VTVM 
and  the  wave  analyzers  simultaneously. 

« 

The  calibration  of  the  accelerometer  was  accomplished  by 
direct  measurement  of  the  piston  face  displacement  at  100  Hz.  and 
at  various  power  settings  with  a  stroboscopic  light  and  a  traveling 
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where  V  is  the  RMS  voltage  reading  on  the  microphone  system  VTVM. 

The  strength  parameter  can  now  be  reformulated  as 

SP  =  $2  VbQ^SjDc2), 

all  of  which  are  observable  or  calculable  quantities.  The  calcu¬ 
lation  or  measurement  of  each  of  these  quantities  introduces  error 
into  the  determination  of  the  strength  parameter;  this  error  is 
estimated  to  be  about  1.0  to  1.5  percent  for  both  the  modes 
investigated. 

D.  FREQUENCY  PARAMETER 

The  frequency  parameter  is  a  quantity  which  indicates  the 
position  of  the  frequency  being  analysed  relative  to  the  fundamental 
resonance  of  the  system.  The  frequency  parameter  is  defined  by 

FP  =  2((0-Wlr)Q1/a)lr 

and  can  be  shown  to  be  related  to  the  experimental  frequency 
increment  by 

FP  =  2(0M0lr)/(C01+-C01J. 

E.  SYSTEM  ALIGNMENT 

The  lower  limit  to  which  significant  measurements  could  be 
taken  was  determined  by  the  distortion  introduced  by  the  piston 
motion.  In  an  effort  to  obtain  meaningful  data  at  the  higher 
harmonics  and  for  any  harmonics  at  low  strengths,  great  care  was 
taken  to  reduce  the  distortion  introduced  into  the  cavity.  The 
distortion  of  the  piston  motion  when  free  of  the  cavity  was  below 
the  general  electronic  noise  level  of  the  accelerometer  system  and 
thus  was  insignificant.  The  bulk  of  the  distortion  came  from 


misalignment  of  the  piston  and  cavity  .orifice.  The  procedure  used 
to  obtain  maximum  alignment  of  the  piston  and  thus  minimize  the 
distortion  was  as  follows:  The  output  signal  of  the  accelerometer 
system,  which  measured  the  piston  motion,  was  displayed  on  the 
dual  trace  oscilloscope  and  this  waveform  was  compared  to  a  pure 
sine  wave.  The  piston  and  exciter  apparatus,  which  was  mounted  in 
an  adjustable  framework,  was  positioned  relative  to  the  cavity 
orifice  until  the  accelerometer  waveform  was  as  clean  and  undis¬ 
torted  as  possible.  In  order  to  minimize  any  effects  which  may 
have  arisen  from  vibration  coupling  of  the  cavity  and  exciter  or 
cavity  and  table,  the  cavity  was  mechanically  isolated  from  both 
the  exciter  anc  the  table.  This  was  accomplished  by  mounting  both 
the  cavity  and  the  piston-exciter  apparatus  on  platforms  con¬ 
structed  of  alternate  slabs  of  wood  and  rubber.  This  restricted 
the  interaction  of  the  cavity  with  the  piston  to  the  piston  ori¬ 
fice  only.  With  proper  alignment,  it  was  possible  to  reproduce 
conditions  where  the  distortion  in  the  input  waveform  to  the 
cavity  was  less  than  0.2%  second  harmonic  and  0.3%  for  third  har¬ 
monic,  as  measured  by  the  H-P  302a  wave  analyzer.  Distortion 
values  for  the  higher  harmonics  were  below  the  general  ambient 
noise  of  the  total  apparatus,  and  thus  were  ignored. 


IV.  PROCEDURES 


A.  PRERUN  PROCEDURES 

Before  any  meaningful  data  could  be  taken,  it  was  important  to 
determine  the  amount  of  vibration  transferred  through  the  cavity 
walls  and  picked  up  by  the  microphone.  To  make  this  determination, 
the  top  of  the  cavity  was  removed,  thus  relieving  the  pressure  seal, 
and  a  wood  plate  added  between  the  piston  and  microphone  to  isolate 
the  microphone  from  the  pressure  waves  induced  in  the  air  by  the 
piston.  As  described  in  the  System  Alignment  section,  the  cavity 
and  the  exciter -piston  assemblies  were  isolated  by  a  floating 
platform  technique.  The  vibration  picked  up  by  the  microphone  was 
measured  at  frequencies  of  interest  and  produced  microphone  system 
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responses  of  8 *10  to  2‘10  volts.  In  all  cases  these  values 
amounted  to  distortion  well  below  0.3%,  which  is  the  maximum 
amount  of  piston  motion  distortion. 

It  was  noted  by  Beech  [ 3 J  that  temperature  changes  produced 
significant  changes  in  the  resonant  frequency  of  the  tube  he  in¬ 
vestigated.  The  temperature  effect  on  the  cavity  investigated 
here  was  very  pronounced.  This  presented  many  problems  as  the 
frequency  of  fundamental  resonance  was  the  reference  for  all  de¬ 
terminations  of  frequency  parameter.  Therefore  it  was  necessary 
to  know  the  value  of  the  fundamental  resonance  at  all  times.  To 
minimize  this  shift  in  the  resonance  frequency,  the  entire  ap¬ 
paratus  was  allowed  to  warm  up  for  a  period  of  one  to  four  hours, 
depending  on  the  strength  parameter  to  be  investigated,  prior  to 


taking  data.  This  procedure  did  not  totally  relieve  the  effect  of 
the  frequency  shift  so  another  calibration  method  was  employed. 

Many  determinations  of  the  fundamental  frequency  were  made  during 
the  taking  of  data,  thus  producing  a  calibration  curve,  see 
Figure  3»  by  which  the  frequency  parameter  values  could  be 
corrected. 

Immediately  prior  to  the  taking  of  each  data  set,  the  align¬ 
ment  of  the  piston  and  cavity  was  checked  as  outlined  in  the  system 
alignment  section  of  Experimental  Considerations.  From  this  point 
onward,  the  procedure  for  data  collection  \vas  dictated  by  the 
quantities  being  measured.  There  was  one  experimental  method  to 
measure  the  Q's  and  E's  of  the  system  in  the  linear  region,  and 
another  to  measure  the  harmonic  content  of  the  non-linear  finite 
amplitude  waves. 

B.  MEASUREMENT  OF  Q(N)  AND  E(N)  IN  THE  LINEAR  REGION 

Data  collection  in  this  region  required  very  little  power, 
but  required  very  fine  control  of  the  power  setting.  This  control 
was  obtained  by  inserting  a  voltage  dividing  circuit  between  the 
frequency  synthesizer  and  the  power  amplifiers.  This  voltage  di¬ 
vider  reduced  the  available  power  by  a  factor  of  ten  and  increased 
the  sensitivity  of  the  power  control  by  the  same  factor  of  ten. 

With  this  arrangement,  a  satisfactory  fine  tuning  paver  control 
was  obtained  at  low  strengths. 

With  the  power  set  at  some  nominal  low  value,  the  frequency 
region  about  the  resonant  or  harmonic  pressure  peak  of  interest 
was  scanned  until  a  maximum  in  pressure  was  indicated  by  the 


microphone  system  response.  The  input  power,  determined  in  the 
linear,  low  strength,  region,  and  kept  constant  during  the  fol¬ 
lowing  measurements.  The  frequency  of  the  driving  piston  was 
changed  above  and  below  the  frequency  of  maximum  pressure,  until 
the  pressure  was  -3  dB  of  the  maximum  pressure  value.  This  deter¬ 
mination  of  the  -3  dB  points  was  made  with  the  302a  wave  analyzer, 
which  measured  the  pressure  amplitude  of  the  cavity  by  way  of  the 
microphone  system. 

It  was  found  that  the  first  seven  overtones  of  the  nOO  mode 
and  the  first  five  overtones  for  the  OoO  mode  were  undisturbed 
by  overlapping  Q  curves.  The  value  of  Q  for  each  response  curve 
was  determined  from  the  frequency  difference  of  the  -3  dB  points, 
divided  by  the  center  point  frequency  between  the  two  -3  dB 
points.  This  center  point  was  assumed  to  be  a  good  estimate  of 
the  pe:'.k  pressure  frequency,  o>  ■ ,  for  the  response  curve  in 
question.  It  was  with  these  values  for  the  system  resonance  fre¬ 
quencies  that  the  E(n)’s  were  determined. 

C.  MEASUREMENT  OF  THE  HARMONIC  CONTENT  OF  FINITE  AMPLITUDE 

STANDING  WAVES 

In  order  to  obtain  finite  amplitude  waves  in  the  cavity,  it 
was  necessary  to  go  to  much  higher  strengths  than  were  required 
for  the  measurement  of  the  Q’s  and  E’s.  With  the  driving  fre¬ 
quency  set  to  some  value  near  the  fundamental,  and  the  peak  pres¬ 
sure  component  of  the  fundamental  held  constant,  the  strength  of 
each  harmon’C  was  measured  by  one  of  the  wave  analyzers.  The 
other  wave  analyzer  monitored  the  fundamental  pressure  ,  which 


v/as  kept  constant  by  varying  the  power  setting  of  the  synthesizer 
as  necessary.  Values  of  the  frequency  parameter  were  chosen  such 
that  the  distortion  curves,  Figures  4  through  9,  could  be  deter¬ 
mined  in  the  areas  of  maximum  distortion.  The  order  in  which 
measurements  were  taken  for  each  distortion  curve  was  randomized 
to  avoid  systematic  error.  Strength  parameters  of  about  0.5,  0.75, 
and  1.0  were  investigated  for  each  mode.  Strengths  below  0.5  did 
not  produce  high  enough  distortion  to  be  measured  with  the  equip¬ 
ment  available.  Strengths  above  1.0  were  not  accessable  due  to 
limitations  in  power . 


V.  RESULTS 


Tables  I  and  II  are  the  data  for  the  linear,  infinitesimal 
amplitude  region.  The  information  on  the  empirical  losses  is  con- 

i 

tained  in  the  Q’s  and  E's.  These  are  the  values  used  in  the  com¬ 
puter  program  to  predict  the  harmonic  distortion  curves  in  Figures 
4  through  9.  The  theoretical  curves  in  these  figures  are  plotted 
along  with  the  experimentally  measured  values  of  harmonic  dis¬ 
tortion.  The  theoretical  predicted  values  were  generated  for  fre¬ 
quency  parameter  intervals  of  0.2  and  were  plotted  as  solid  curves. 
The  experimentally  measured  values  were  plotted  au  dots  where  they 
occured.  Data  were  taken,  and  theoretical  predictions  made,  for 
strength  parameters  of  0.513>  0.782,  and  1.00  for  the  nOO  mode  and 
0.507,  0.780,  and  1.00  for  the  OnO  mode. 

The  point  of  intersection  of  the  error  bars  on  the  P  /P1  theo- 

retical  curve  in  Figure  4  corresponds  to  the  value  of  EP  =  2Q  E(2) 

0)  1 
2r 

where  to  =  .  It  is  important  to  note  that  the  n  =  2  distortion 

curves  peak  when  the  system  is  driven  at  this  frequency.  That  is, 
when  the  driving  frequency  to  is  equal  to  the  classical  fundamental 
of  C02r,  there  is  maximum  content  of  P2<  The  point  where  this  oc¬ 
curs  for  each  P2/P^  curve  is  indicated  by  arrow?  in  Figures  5 
through  9  and  by  crossed  error  bars  in  Figure  4.  The  error  in  the 
determination  of  these  points  comes  mainly  from  the  error  in  E(2) 
that  appears  in  the  above  equation  for  FP.  The  error  in  is  of 
much  less  significance  as  is  indicated  in  Table  III.  The  error  in 
the  P2/P^curves  in  Figures  5  through  9  are  approximately  the  same 
as  those  indicated  in  Figure  4» 


Table  1 


nOO  Mode 


n 

f  JH2) 

f+(H2) 

Af(Hz) 

VH2> 

Q 

E  x  103 

1 

677.80 

679.70 

1.90 

678.75 

357.24 

0 

2 

1365.74 

1368.49 

2.75 

1367.12 

497.13 

7.0 

3 

2059.17 

2062.84 

3.67 

2061.01 

561.58 

12.0 

4 

2787.-37 

2792.73 

5.36 

2790.05 

520.53 

27.5 

5 

3423.00 

3427.80 

4.80 

3425.40 

713.63 

9.1 

6 

4126.41 

4131.91 

5.50 

4129.16 

750.76 

13.6 

7 

4803.92 

4809.93 

6.01 

4806.93 

799.82 

11.4 

l 

678.02 

679.95 

1.93 

678.99 

351.81 

0 

OnO 

n 

Mode 

f_(Hz) 

f+(H2) 

Table  II 

Af (Hz)  fn(Hz) 

Q 

E  x  103 

1 

813.99 

816.08 

2.09 

815.04 

9. 97 

0 

2 

1633.16 

1636.39 

3.23 

1634.78 

506 . 12 

2.8 

3 

24l5.ll 

2420.57 

5.46 

2417.84 

442, 

"11.3 

4 

3143.74 

3150.10 

6.36 

3146.92 

494.80 

-34.9 

5 

3911.05 

3920.35 

9.30 

3915.70 

423.62 

-39.4 

l 

814.23 

816.34 

2.11 

S15.2Q 

386.39 

0 

25 


incy  Parameter 


Table  III  contains  the  results  of  five  separate  determinations 
of  the  Q’s  and  E’s  taken  in  the  infinitesimal  amplitude  region. 
These  measurements  were  made  to  provide  an  estimate  of  the  error 
in  the  theoretical  curves.  In  all  cases  the  values  used  to  gene¬ 
rate  the  curves  were  very  close  to  the  average  values  determined 
later  and  the  regeneration  of  the  curves  was  therefore  unwarranted. 
Table  III  contains  both  sets  of  values  for  direct  comparison. 
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Table  III 


nOO  Mode 


Values  Used  to 

Average  Values  Determine  Curves 


n 

Q(n) 

E(n) ‘103 

Q(n) 

E(n) • 103 

1 

358.2  +  0.4 

357.2 

2 

498.5  +  4.7 

6.8  +  o.l 

497.1 

7.0 

3 

559.8  +  0.8 

11.6  +  0.2 

561.6 

12.0 

OnO  Mode 


l 

359-8+5-3  . 

390.0 

— 

2 

515.7  +  5.1  2.6  +  0.1 

509.1 

2.8 

3 

459.2  +  7.5  -11.3  +  0.1 

442.8 

-11.3 
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VI .  CONCLUSIONS 


It  appears  that  the  present  theoretical  model  successfully 
predicts  the  major  features  of  the  harmonic  content  for  finite- 
amplitude  standing  waves  in  a  rigid-walled  cavity.  In  all  cases 
investigated,  the  shape  and  general  magnitude  were  predicted;  see 
Figures  4  through  9. 

It  is  interesting  to  note  the  position  of  each  theoretical 
curve  relative  to  its  corresponding  experimental  curve.  For  the 
nOO  mode,  the  theoretically  predicted  values  are  consistently 
above  the  measured  values.  The  situation  for  the  OnO  mode  is  the 
same  for  P  /P  but  opposite  for  P  /P  .  A  possible  source  of  this 

O  1  2  i 

difference  is  the  large  difference  in  the  geometry  encountered  by 
each  mode.  The  nOO  mode  is  a  standing  wave  whose  wave  vector  is 
perpendicular  to  the  wall  with  the  microphone  orifice,  which  must 
be  considered  a  very  small  perturbation  due  to  its  small  size.  On 
the  other  hand,  the  OnO  mode  is  a  standing  wave  with  its  wave 
vector  perpendicular  to  the  piston  orifice,  which  is  a  large  per¬ 
turbation  in  the  geometry.  In  addition,  the  fact  that  the  wave 
generating  surface,  the  piston  face,  is  normal  to  the  nOO  wave, 
vector  and  perpendicular  to  the  OnO  vector  could  cause  a  con- 
sistant  difference  in  the  behavior  of  the  modes.  It  is  certainly 
unclear  to  this  investigator  just  what  the  specifics  are  that 
explain  the  trends  in  the  theory/exper iraent  relationship.  It  is 
clear  however,  that  the  theory  of  Coppens  and  Sanders  as  modified 
to  include  empirical  losses,  does  predict  the  essence  of  the 
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behavior  of  finite-amplitude  standing  waves  in  the  rigid-walled 
cavity  investigated. 
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APPENDIX  A 


In  the  investigation  by  Winn  [2],  it  was  found  that  the  tubes 
analyzed  did  not  behave  in  a  manner  consistent  with  Rayleigh- 
Kirchoff  wall  losses.  In  an  effort  to  determine  if  the  interior 
finish  of  the  tubes  was  a  critical  factor  that  determined  the  loss 
mechanisms,  the  six  foot  tube  was  reanalyzed  after  refinishing  the 
interior  to  a  highly  polished  state.  These  results  were  compared 
to  the  results  for  the  tube  in  a  smooth  but  not  polished  condition 
as  determined  by  Winn. 

It  was  found  that  the  absorption  coefficients,  defined  as 


CO  -  a) 
n+  n 


were  not  significantly  changed  from  the  values  found  by  Winn,  and 
both  sets  of  values  obeyed  the  Rayleigh-Kirchoff  predictions 
indicated  by  the  solid  line  in  Figure  10.  The  behavior  of  the 
overtones  (for  the  polished  tube)  was  found  to  be  significantly 
different  (from  the  overtone  behavior  in  the  unpolished  tube), 
and  was  significantly  different  from  the  Rayleigh-Kirchoff  pre¬ 
dictions  as  indicated  by  a  solid  line  in  Figure  11. 


BIBLIOGRAPHY 


1.  Coppens,  A.B. ,  and  Sanders,  J.V.,  "Finite-amplitude  Standing 
Waves  in  Rigid  Walled  Tubes,"  J.  Acoust.  Soc.  Am.,  V.  43, 

P.  516-529,  March  1968. 

2.  Winn,  J.R.,  Fourier  Analysis  of  Experimental  Finite-amplitude 
Standing  Waves ,  Thesis,  Naval  Postgraduate  School,  Monterey, 
California,  1971. 

3.  Beech,  W.L.,  Finite-amplitude  Standing  Waves  in  Rigid-walled 
Tubes,  Thesis,  Naval  Postgraduate  School,  Monterey,  California, 

19^77 

4.  Kirchoff,  G.,  Ann.  Phvs.  Leipzig,  V.  134,  P.  177-193,  1868. 

5.  Lamb,  L.  ,  Dynamical  Theory  of  Sound,  2d  ed. ,  Chap.  VI, 

Edward  Arnold  and  Co.,  London  England,  1925. 

6.  Fay,  R.D. ,  "Successful  Method  of  Attack  on  Progressive  Finite 
Waves,"  J.  Acoust.  Soc.  Am.,  V.  28,  P.  910-914,  September 
1956. 

7.  Blackstock,  D.T.,  "Propagation  of  Plane  Sound  Waves  of  Finite 
Amplitude  in  Nondissipative  Fluids,"  J.  Acoust .  Soc.  Am. , 

V.  34,  P.  9-30,  January  1962. 

8.  Lord  Rayleigh,  Theory  of  Sound,  2d  ed.,  V.  1  and  2,  Dover 
Publications,  Inc.,  New  York,  1945. 

9.  Fay,  R.D.,  "Plane  Sound  Waves  of  Finite  Amplitude,"  J.  Acoust. 
Soc.  Am..  V.  j,  P.  222-241,  1931. 

10.  Blackstock,  D.T.,  "Convergence  of  the  Keck-Beyer  Perturbation 
Solution  for  Plane  'Waves  of  Finite  Amplitude  in  a  Viscous 
Fluid,"  J.  Acoust.  Soc.  Am..  V.  39,  P.  411-413,  1966. 

11.  Fox,  F.E.,  ana  Wallace,  W.A.,  "Absorption  of  Finite  Amplitude 
Sound  Waves,"  J.  Acoust.  Soc.  Am..  V.  26,  P.  994-1006,  1956. 

12.  Gol :berg,  Z.A.,  Akust.  Zhur .  V.  3,  P.  149,  1957,  (English 
Trans.;  Sov.  Phys .  -  Acoust.  V.  3,  P.  157,  1957). 

13.  Keck,  W.,  and  Beyer,  R.T.,  "Frequency  Spectrum  of  Finite 
Amplitude  Ultrasonic  Waves  in  Liquids,"  Phys.  Fluids,  V.  3, 

P.  346-352,  I960. 

14.  Keller,  J.B.,  "Finite  Amplitude  Sound  Produced  by  a  Piston 
in  a  Closed  Tube,"  J.  Acoust.  Soc.  Am.,  V.  26,  P.  253-254, 

1954. 


39 


15.  Weston,  D.E.,  Proc.  Phys .  Soc.,  (London)  B66,  695-709,  1953. 

16.  Ruff,  P.G.,  Finite  Amplitude  Standing  Waves  in  Rigid  Walled 
Cavities ,  Thesis,  Naval  Postgraduate  School,  Monterey, 
California,  1967. 

17.  Coppens,  A.B.,  "Theoretical  Study  of  Finite-Amplitude 
Traveling  Waves  in  Rigid-Walled  Ducts:  Behavior  for  Strengths 
Precluding  Shock  Formation."  J.  Acoust.  Soc.  Am.,  V.  49, 

P.  306-318,  1969. 

18.  Coppens,  A.B.,  personal  communication. 


» 


40 


